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ABSTRACT 

Our discovery of ISWASP J093010.78-I-533859.5 as a probable doubly eclipsing quadruple system, containing a contact binary with 
P ~ 0.23 d and a detached binary with P ~ 1.31 d, was announced in 2013. Subsequently, Koo et al. confirmed the detached 
binary spectroscopically, and identified a fifth set of static spectral lines at its location, corresponding to an additional non-eclipsing 
component of the system. Here we present new spectroscopic and photometric observations, allowing confirmation of the contact 
binary and improved modelling of all four eclipsing components. The detached binary is found to contain components of masses 
0.837 ± 0.008 and 0.674 ± 0.007 Mq, with radii of 0.832 ± 0.018 and 0.669 ± 0.018 Rq and effective temperatures of 5185;^2o 
4325+J5 K, respectively; the contact system has masses 0.86 ± 0.02 and 0.341 ± 0.011 Mq, radii of 0.79 ± 0.04 and 0.52 ± 0.05 Rq, 
respectively, and a common effective temperature of 4700 ± 50 K. The fifth star is of similar temperature and spectral type to the 
primaries in the two binaries. Long-term photometric observations indicate the presence of a spot on one component of the detached 
binary, moving at an apparent rate of approximately one rotation every two years. Both binaries have consistent system velocities 
around -11 to -12 km s ', which match the average radial velocity of the fifth star; consistent distance estimates for both subsystems 
of d = 78 ± 3 and <7 = 73 ± 4 pc are also found, and, with some further assumptions, of <7 = 83 ± 9 pc for the fifth star. These findings 
strongly support the claim that both binaries - and very probably all five stars - are gravitationally bound in a single system. The 
consistent angles of inclination found for the two binaries (88.2 ± 0.3°and 86 ± 4°) may also indicate that they originally formed by 
fragmentation (around 9-10 Gyr ago) from a single protostellar disk, and subsequently remained in the same orbital plane. 
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1. Introduction 


During the course of a search for orbital period variations 
in short- period eclipsing bina ry candidates i n the Supe r WASP 


[- period eclipsing bin^ ^ 

arc hive (j Pollacco et alJ l20()^ . described in iLohr et al.l d2012ll 
and lLohr et al. ( 2013bl) . several unusual systems of particular as¬ 
tronomical interest were encounter ed. Amongst these , IS WASP 
J234401.81-212229.1, discussed in lLohr et al](l2013ah and iKoenI 
(12014 . appears to be a triple system containing an M-hM 
dwarf contact binary. Here we explore ISWASP J093010.78- 
-1-533859.5 (hereafter J093010), an apparent quintuple system 
containing two eclipsing binaries, one of only a handful of 
known doubly eclipsing multiples. 


Higher-order multiple systems are of value in general for 
testing models of stellar formation and long-term dynamical sta¬ 
bility. When they contain two double-lined spectroscopic and 
eclipsing binaries, they provide a highly unusual opportunity to 
determine the component stars’ physical parameters and so to 
understand the structure and potential origins of the whole sys¬ 
tem. 


2. Background 


A near neighbour of J093010 (ISWASP J093012.84-H533859.6, 
hereafter J093012) was identified as a candidate short-period 
eclipsing bin ary in iNorton et al.l (1201 ll) . and its orbital period 
confirmed in iLohr et aT~(l2012l) . However, in our subsequent 
more thorough search for eclipsing candidates, J093010 was 
identified as exhibiting the same period and light curve shape, 
but a higher mean flux a nd greater amplitud e of flux variability; 
hence, it was regarded in lLohr et al.l (1201 3bh as the probable true 
source of the observed eclipsing variation, with J093012 being a 
fainter near-duplicate light curve, falling within the same photo¬ 
metric aperture used by SuperWASP. At the location of J093012 
there is a magnitude 18 source listed in the USNO-Bl catalogue, 
separated from J093010 by approximately 18", which would not 
have been detectable by SuperWASP in its own right. 


The eclipsing source’s apparently significant period change, 
associated with an erratic O-C diagram, had been r ejected ini¬ 
tially as the r esult of con t aminati on by a nearby star (iLohr et al.l 
I20l4 . but in iLohr et all (l2013bt) a fuller explanation was pur¬ 
sued. Prior to analysis, the data for J093010 and J093012 were 
combined to maximize the available observations. J093010’s 
light curve, folded at 19674.574 s, then showed a typical con- 
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tact binary shape, but with significant non-Gaussian data scatter 
below the main curve. A visual examination of the object’s full 
light curve suggested the cause was additional deep eclipses on 
certain individual nights, implying a second eclipsing body in 
the field of view. A frequency power spectrum also supported 
an additional periodic signal near 1.3 d. Subtracting the median 
binned light curve (corresponding to the contact binary) from the 
data revealed the light curve of an EA-type (detached) eclipsing 
binary with period 112799.109 s. 

The juxtaposition of these two binaries on the sky did not 
seem to be coincidental. Two sources had been observed at this 
location by Hipparcos as TYC 3807-759-1 and TYC 3807-759- 
2, with equivalent Johnson V magnitudes 9.851 and 10.990, re¬ 
spectively; a separation of T.'88; and a proper motion of pmRA: 
-8.0 mas yr', pmDE; -9.4 mas yr', measured at their joint pho¬ 
tocentre. The tiny angular separation and similar magnitudes 
favoured a plausible interpretation of the two eclipsing binary 
systems as being gravitationally bound in a quadruple doubly 
eclipsing system. 

Only five other doubly eclipsing quadruple systems had 



J093010 had a shorter period than any system in the other five 
quadruples, making it particularly amenable to further observa¬ 
tions. 

iKoo et all ( 1201 4|) obtained BV photometry over several 
months in 2012-2013 for each eclipsing binary separately, and 
time-series spectra for the whole system. Their light curves made 
clear that the brighter signal (TYC 3807-759-1) was associated 
with the detached binary, and they consequently termed this 
source J093010A, and the fainter source J093010B (containing 
the contac t binary). We follow this convention here. Although 
iKoo et all were unable to detect a spectral signature of the as¬ 
sumed contact binary in their combined spectra (due in part to 
the length of their exposures), they did observe clear line split¬ 
ting and shifting with the period of the proposed detached binary, 
confirming its reality. They also unexpectedly detected an addi¬ 
tional set of static spectral lines which they interpreted as a third 
component in J093010A, which would make the whole system 
a highly unusual doubly eclipsing quintuple. 

Independently, we obtained time-series spectroscopy over 
three nights in 2012-2013 for J093010A and J093010B sepa¬ 
rately, and near-simultaneous RGB photometry for the whole 
system to help establish the phases of the obs ervations. Our new 
data sets thus complement those provided in iKoo et all (l2014l) . 
allowing us to model both binaries quite fully; the revisited 
SuperWASP observations also provide further long baseline in¬ 
formation. 

2.1. Spectroscopy 

49 long-slit spectra were obtained by E. Gillen for J093010 on 
22 and 31 December 2012, and on 2 January 2013, using the red 
arm of the ISIS spectrograph on the 4.2 m William Herschel 
Telescope at La Palma. The R1200R grating provided an in¬ 
termediate resolution of 0.24 A per pixel. A usable wavelength 
range of ~8350-9000 A was obtained, to capture the Ca II triplet 
important in low-mass stars. Exposure lengths ranged from 5 to 
360 s (see Table [1]): short enough in all cases to avoid motion 
blur for the contact binary. 


Auerage of lines 2301 to 3300 of rl856221.fits 
ISUASP 



Fig. 1. Left: Section of spectra for J093010A and J093010B 
from a single 360 s exposure, showing spatial separation. Right: 
Plot of spectral profile from same exposure, cut perpendicular to 
the dispersion axis, showing two clear peaks in the line disper¬ 
sion function. The stronger peak on the right-hand side corre¬ 
sponds to J093010A. 

The images were flat-fielded and bias-corrected by E. Gillen 
as part of a larger observing programme. On all but one image, 
two partially overlapping dispersion lines were clearly visible 
(e.g. Eig. [TJ, corresponding to J093010A and J093010B, and 
we were able to extract separate spectra for each, using stan¬ 
dard IRAE tools and optimal extraction. Wavelength calibration 
was carried out using arc spectra taken with CuArNe lamps. S/N 
ratios reached 340 for the brighter J093010A and 220 for the 
fainter J093010B in the longest exposures, when both systems 
were near quadrature, and fell to 45 and 25 respectively in the 
initial 5 s exposure. 

Preliminary radial velocities were extracted using the IRAE 
task EXCOR, using a synthetic spectrum with T =4500 K and 
log ^=4.5 as the comparison template. These allowed approx¬ 
imate determination of mass ratios and radial velocity curve 
amplitudes, which provided initial input parameters fo r spectral 
disen tangling of J093010A and B. The KOREL code (iHadraval 
1201 2h . implemented on the Virtual ObservatorjQ, w as used t o 
disentangle sections of the spectra prepared with iHadraval ’s 
PREKOR software, giving improved values for mass ratio and 
radial velocities. The disentangled spectra were also compared 
with appropriately broadened template spectra extracted from 
the PJJOENIX synthetic stellar library (llJusser et al.ll2013lR to 
estimate the temperatures of individual components. 

2.2. Photometry 

67 images of J093010 (55 each in Baader G and B filters, 57 
in R) were taken with the Open University’s robo tic 0.425 m 
PIRATE telescope (iHolmes et al.ll20lll lKolbll2014ll in Mallorca 
on the night of 30-31 December 2012, monitored by R. Busuttil. 
All exposures were 80 s; short enough to avoid phase smearing 
for the shorter-period system. Corrections for bias level, dark 
curi'ent and flat-fielding were made to the frames using standard 
IRAE tools, and four comparison stars were identified on the 
frames: TYC 3807-1509-1, TYC 3807-54-1, TYC 3807-1503- 
1 and TYC 3807-621-1. These had catalogue colours indicative 
of classes between late E and G, and V magnitudes between 9 
and 10, similar to J093010. They were checked for short and 
long-term variability with their SuperWASP light curves, and did 

* https://stelweb.asu.cas.cz/vo-korel 

^ phoenix.astro.physik.uni-goettingen.de 
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not exhibit significant variations. Aperture photometry was car¬ 
ried out on all the stars using the IRAF APPHOT package, and 
the light curves of the comparison stars were combined in IDL. 
Differential light curves were then obtained for J093010 relative 
to this combined comparison curve. 

The Super WASP archive data for J093010 were also recon¬ 
sidered; there are 5964 photometric points in V, observed be¬ 
tween October 2007 and April 2009. (The fainter duplicate ob¬ 
ject J093012 has 5950 points covering the same time span; its 
data were used here as a check on the results found for J093010.) 
For reference, SuperWASP exposure lengths are 30 s. Initial val¬ 
ues for the orbital periods of the two ecli psing binaries wer e 
found with a custom IDL code described in lLohr et alJ (l2014bh : 
19674.594+0.005 s and 112799.10+0.15 s; then a second code 
was written to separate the two eclipsing signals. The light curve 
was first folded on the shorter period corresponding to the con¬ 
tact binary, and phase-binned to give a smooth mean curve (90 
bins were used); an optimally-weighted average was then found 
for the data points in each bin (which corresponded closely to the 
visible contact binary signal). A spline curve was interpolated to 
these binned average points, and subtracted from the full light 
curve. The residue was then folded on the longer period cor¬ 
responding to the detached binary, and a binned average curve 
obtained as before. This was again subtracted from the full light 
curve, to leave a clean contact binary light curve as output. This 
clean mean contact curve was subtracted a second time from the 
full data set, to leave a clean detached binary light curve for out¬ 
put. Our initial period-searching code was run again on the sep¬ 
arated files, and its resulting period values were used to initialize 
the signal-separation code for a second iteration. This continued 
until convergence was reached for both periods. 

The BV light curves of iKoo et al.l (1201 4|) were also used in 
our modelling; most details of the photometry are given in that 
paper, though we add that the exposure lengths were below 25 s, 
avoiding any risk of phase-smearing for even the shorter-period 
binary. 

3. Results 

3.1. Orbital periods 

The orbital periods based on SuperWASP data rapidly converged 
on 19674.47+0.03 s for the contact binary and 112798.90+0.16 s 
for the detached binary i.e. 0.2277138(3) d and 1.3055428(19) d 
respectively. (The periods for J093012 did not converge, but 
oscillated within the narrow ranges 19674.43-19674.52 and 
112798.9-112800.0, which include and support the results for 
J093010. ) These pe riods fall within the uncertainties of those 
found bv lKoo et al.l from their photometry; 0.2277135(16)d and 
1.30550(4) d, though their figure for the detached binary’s pe¬ 
riod is based on a single measured primary eclipse time and three 
secondary eclipses. 

No period change was detected within the SuperWASP data 
for either binary; extrapolating forward from the most precisely 
measured times of minimum for the SuperWASP data, using our 
periods, gives an O-C value of just -30 s for the lKoo et aLm ri- 
mary minimum for the detached system, and —1-200 s for their 
primary minima for the contact system. The small size of these 
discrepancies further supports both the absence of significant pe¬ 
riod change (especially in the detached system) and the reliabil¬ 
ity of our periods, calculated over long base lines using whole 
light curves, rather than minimum timings alone. We therefore 
prefer our periods for the remainder of the analysis, but use 



Fig. 2. Extract of WHT spectra for J093010A at approximate 
phases 0.16, 0.64 and 0.75 (top to bottom). Relatively isolated 
absorption lines at around 8793, 8806 and 8823 A exhibit three- 
way splitting, with the outer pair also shifting between phases. 
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Pixel Shift 

Fig. 3. Cross-correlation plot for J093010A at phase 0.75, ex¬ 
hibiting three strong peaks corresponding to three stellar com¬ 
ponents (primary on the right, secondary on the left, tertiary in 
the middle). The broadest absorption lines were excluded from 
the cross-correlation to minimize the radial velocity uncertain¬ 
ties. 

iKoo et all ’s more precise and recent primary minimum timings 
to calculate phases. 

3.2. Radial velocities and spectroscopic parameters 

IKoo et all measured radial velocities for J093010A by fitting 
three Gaussians simultaneously to five isolated absorption lines 
in their spectra, in the range 6400-6800 A. Our longer wave¬ 
length range contained such a profusion of blended lines that this 
approach would have been impossible; indeed, it was hard to de¬ 
tect three-way line splitting from visual examination of the spec¬ 
tra over most of the range (Fig.|2]illustrates a small region where 
the triplets are relatively isolated). However, cross-correlation 
revealed three extremely clear correlation peaks for every spec- 
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Phose 

Fig. 4. Radial velocity curves for J093010A. The primary’s mea¬ 
surements are indicated by squares, the secondary by diamonds 
and th e third component by triangles. Data from iKoo et afl 
(l2014ll is in fainter print (orange and light blue in the online ver¬ 
sion of the paper) and our new results ar e in bold print (red, 
dark blue and green). Phase-folding uses iKoo et al.l ’s primary 
minimum time of HJD 2456346.78443 and our optimum period. 
PHOEBE model fits are overplotted (solid black curves) and the 
location of the modelled system velocity is shown by the dashed 
line. 

trum (e.g. Eig. [3, allowing preliminary radial velocities to be 
measured easily for the system. 

Spectral disentangling was then achieved for three subre¬ 
gions of our spectra containing narrow and well-defined lines, 
and converged on consistent values for mass ratio {q - 0.806 + 
0.007) and semi-amplitude of the primary’s radial velocity curve 
{Ki = 99.4 + 0.7 km s '). The resulting radial velocities were 
corrected to heliocentric values by disentangling telluric lines 
observed in the region 8250-8320 A, and are given in Table [1] 
(These are very close to those found by the cross-correlation 
method, but with sma ller uncerta inties.) Eig.|4]shows our results, 
together with those of iKoo et~al] : where our observations overlap 
with theirs (around phases 0.15 and 0.65) there is close agree¬ 
ment in both the amplitudes and absolute values of the curves; 
we have also fortuitously been able to obser ve the syste m around 
its secondary maximum (phase 0.75) where l^o et al.l had a gap. 

Also notable is our confirmation of the third stro ng and near - 
static component of the detached system spectra; iKoo et alf s 
fifth star. Since we have obtained separate spectra for J093010A 
and J093010B, we can be certain that this additional source oc¬ 
curs at the same location on the sky as the detached eclipsing 
binary, rather than being near the contact system or mid-way be¬ 
tween the two. There is also very good reason to regard it as 
part of a gravitationally bound system including the detached bi¬ 
nary; its radial velocities (average -11.6 + 1.5 km s ') are visibly 
very close to the cross-over system velocity of the binary (mod¬ 
elled as -12.3 + 0.2 km s '). While longer-term radial velocity 
and astrometry measurements would be required to demonstrate 
convincingly that the fifth star forms a triple subsystem with the 
detached binary (rather than, for example, the two binaries being 
more closely bound to each other), the coincidence of its spec¬ 
trum with that of the binary in J093010A (see Eig.[TJ makes this 
a more plausible hierarchy. 

We were also able to measure radial velocities from most 
of the spectra for the candidate contact system in J093010B. 
Here, splitting and shifting of the strongest lines (primarily the 



Wavelength (angstroms) 

Fig. 5. Extract of WHT spectra for J093010B at approximate 
phases 0.24, 0.66,0.88 and 0.99 (top to bottom). The Call triplet 
of absorption lines at around 8498, 8542 and 8662 A exhibit 
clear splitting and shifting between phases. 
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Fig. 6. Cross-coiTelation plot for J093010B at phase 0.24, ex¬ 
hibiting two strong peaks corresponding to the binary compo¬ 
nents (primary on the left, secondary on the right). The two 
broadest calcium lines at 8542 and 8662 A were excluded from 
the cross-correlation to minimize the uncertainty in the radial 
velocity measurements. 


Call triplet) was apparent from visual inspection of the spec¬ 
tra alone (e.g. Fig. |5]l, confirming this system as a double-lined 
spectroscopic binary also. The results from cross-correlation 
(e.g. Fig. |6]l again provided starting values for spectral disen¬ 
tangling, and here we were able to disentangle the full spec¬ 
trum as a single region, giving much smoother radial velocity 
curves with much smaller uncertainties than had been achieved 
by cross-correlation (Table[T]i. A mass ratio of 0.397 + 0.006 and 
Ki - 105.1 + 1.3 km s ' were converged upon by disentangling; 
a slightly different value of reference minim um was als o implied 
by the spectra; HJD 2456288.87687, where l^o et al.l ’s nearest 
minimum (based on photometry) was HJD 2456288.87879. 

Our results are shown in Fig.|2l we have managed to capture 
the system around its primary maximum (phase 0.25) and on 
each side of its secondary maximum (near phases 0.6 and 0.9), 
giving sufficient phase coverage for decent modelling (though 
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Fig. 7. Radial velocity curves for J093010B, based on our WHT 
observations. The primary is shown with red squares and the sec¬ 
ondary with blue diamonds; uncertainties are generally smaller 
than the symbol size. Phase-folding uses HJD 2456288.8780, 
found by PHOEBE modelling of radial velocity and light curves 
simultaneously, and our optimum period. The dashed line shows 
the location of the (modelled) system velocity. 
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Wavelength / Angstroms 

Fig. 8. Disentangled renormalized spectrum for primary compo¬ 
nent of J093010B (black solid lines), with best-fitting PHOENIX 
synthetic template overplotted (red dotted lines). The prominent 
absorption lines are Ca II, Mg I, Ee I, Ti I and AI I. 

further observations around phases 0.4 and 0.75 would still be 
desirable). The cross-over velocity of the two curves is -11.3 + 
0.7 km s"': very close to the detached binary system velocity and 
the third component velocity in J093010A, strongly supporting 
the hypothesis that the two binaries and the fifth star are gravita¬ 
tionally bound within a quintuple system with a common motion 
relative to the Sun. 

The disentangled spectrum for the primary component of 
J093010B was best fitted by a PHOENIX synthetic template 
with T = 4700 + 50 K and log g=4.5 (Eig. [D. The sec¬ 
ondary component’s spectrum was of very similar shape, but 
had far lower S/N and so was not used to assess the common 
temperature of the contact binary. The temperature found here 
spectrosco pically is entirely consistent with that suggested by 
iKoo et akl ’ s B-V colour for J09301 OB during secondary eclipse 
i.e. 4680 + 50 K. 

The three components of J093010A were difficult to dis¬ 
entangle convincingly owing to heavy blending of strong lines 


Fig. 9. Part of disentangled renormalized spectrum for tertiary 


component of J093010A (see caption to Fig. [8]). 
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Fig. 10. SuperWASP folded light curve for J093010B, converted 
to V magnitudes. The maximum magnitude is that of the quin¬ 
tuple system as a whole. 

and our rather limited phase coverage; however, a reasonably 
good fit to the tertiary component was achieved, since it exhib¬ 
ited narrower, deeper lines than the detached binary components, 
presumably owing to greater rotational broadening in the latter. 
Templates with T = 5100 + 200 K and log g=4.5 (Fig. |9]) pro¬ 
vided the best matches, suggesting an early K spectral type for 
this star. Measurements of the equivalent widths of several dis¬ 
tinct triplets of lines in the original spectra also indicated that 
the primary’s lines were somewhat stronger than those of the 
tertiary, which could imply a slightly higher temperature (and 
mass). Spectroscopy with better phase coverage would be desir¬ 
able in future to allow direct determination of the temperatures 
of all three components in J093010A. However, we were able 
to estimate temperatures for the detached binary components 
from the photometry, during the modelling process described in 
Sect. El 

3.3. Light curve characteristics 

Moving to the ph otometric results, we can compare the light 
curves obtained bv iKoo etaP (l2014t) for J093010A and B sepa¬ 
rately with the SuperWASP light curves extracted from the com¬ 
bined light variation for the whole system. In each case, the 
SuperWASP amplitudes are smaller, and the magnitudes lower. 
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Fig. 11. SuperWASP folded light curve for the detached binary 
in J093010A, converted to V magnitudes. The maximum mag¬ 
nitude is that of the quintuple system as a whole. 


Fig. 13. Normalized differential PIRATE light curves for 
J093010 in Baader R (red, top), G (green, middle) and B (blue, 
bottom) hlters. 
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Fig. 12. iKoo et al.1 light curves for the detached binary in 
J093010A, in V (top) and B (bottom), converted from differen¬ 
tial instrumental magnitudes t o standard magnitudes using nor¬ 
malization indices provided in iKoo et'aP (their Table 6). 

because there is an effective third light included containing the 
(maximum) contributions from the other three stars. For the con¬ 
tact binary the SuperWASP curve (Fig.fTOli is evidently the same 
general shape as the BV curves (Fig. 2 in iKoo et al.h : the sec¬ 
ondary eclipses are visibly flat-bottomed with the primary mini¬ 
mum fractionally deeper than the secondary; the secondary max¬ 
imum is also slightly lower t han that of the primary, though this 
effect is perhaps stronger in iKoo et al.l ’s curves, which would 
support the presence of a long-term static spot on one of the 
components. 

For the detached binary light curves in J093010A (Figs.fTTI 
and[T2]l, the relative dept hs of the tw o eclipses are comparable in 
the SuperWASP data and iKoo et al.l ’s data, and together with the 
sharp ingresses and egresses of the eclipses they support a well- 
detached EA-type system. The eclipses are equally spaced and of 
the same duration, indicating no signihcant orbital eccentricity. 
The primary eclipse appears flat-bottomed while the secondary 
eclipse’s shape is arguably more curved; if taken at face value 
this would suggest that the secondary’s radius is (surprisingly) 
greater than that of the primary. The out-of-eclipse region would 
be expected to be essentially flat in an EA-type light curve, but 
is not. In the SuperWASP curve, the region between phases 0.1 


and 0.4 is far more scattered and noisy than the region between 
phases 0.6 and 0.9; moreover, the maximum is near the primary 
m inimum, w hile it occurs just before the secondary minimum 
in iKoo et'Sl ’s curves; they suggest that this latter feature results 
from a moving spot on one of the components. 

To explore this idea further, we split up the SuperWASP 
detached binary curve by time, into three sections of contigu¬ 
ous data covering several months. The prohle of the out-of- 
eclipse region clearly changed over time, with its minimum oc¬ 
curring around phase 0.8 near HID 2454510, phase 0.4 near 
HID 2454820 and phase 0.3 near HID 2454870 i.e. apparently 
moving negatively in phase at a rate of (very) roughly one cycle 
every two years. The superposition of these three (fairly smooth) 
curves produced the apparent scatter in the curve of Fig.[TT] The 
most plausible explanation for this moving prohle would indeed 
seem to be a spot on the surface of one of the detached system 
components, rotating at a slightly different rate from the star as 
a whole, perhaps because of differential rotation. 

Finally, Fig. [13] shows the PIRATE RGB photometry ob¬ 
tained for the whole system during one night. As predicted by 
the SuperWASP ephemeris, contact binary primary minima oc¬ 
cur around HJD 2456292.519 and 2456292.748, while the two 
eclipsing systems’ secondary minima overlap in between (this is 
most clearly seen in the G curve, where the detached system sec¬ 
ondary minimum occurs around HJD 2456292.60 and the con¬ 
tact one around HJD 2456292.63). The second night of our WHT 
spectra begins just as this photometry ends, and covers phases 
~0.07-0.25 of the contact system, which is entirely consistent 
with the light curve. Additional nights of photometry would be 
useful to establish the average shapes of the three-colour light 
curves; however, we note that the three curves have very simi¬ 
lar shapes. The primary eclipses of J093010B are perhaps frac- 
tionally d eeper in B than R, with G intermediate. Similarly, 
iKoo et all s standardized BV magnitudes for J093010A indicate 
that the primary eclipses are deeper in B, while the secondary 
eclipses are deeper in V i.e. the primary is hotter than the av¬ 
erage of J093010A’s three components, while the secondary is 
cooler. 

3.4. Flux contributions 

We can make further estimates of the relative flux contribu¬ 
tions of the hve stars from the photometry and spectroscopy. 
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Starting with the SuperWASP light curves in SuperWASP flux 
units, if we assume that J093010A and J093010B’s angles of 
inclination are close to 90° (since both systems contain flat- 
bottomed eclipses supporting near-totality), the depth of each 
component’s eclipse gives a minimum value for its flux contribu¬ 
tion (J093010A primary: 38, secondary: 15; J09301 OB primary: 
25, secondary: 23). Adding these four fluxes and subtracting 
them from the total flux (160) when no component is eclipsed 
gives us a maximum flux estimate for the fifth component (59). 
These values would suggest that the tertiary makes up at most 
53% of the flux in J093010A (in V), and t he primary and sec¬ 
ondary at least 34% and 13% respectively. iKoo et aLf s V band 
data similarly supports contributions of around 53%, 37% and 
10% respectively, using the same approach; in B the proportions 
are 53%, 40% and 8%. 

Furthermore, the SuperWASP fluxes allow us to estimate 
J093010B:J093010A flux contributions at the former’s maxi¬ 
mum and primary minimum as 43% and 21%. These figures can 
be compared with the spectroscopic fluxes, found from fits to the 
continua of the extracted spectra, evaluated at 8500 A and scaled 
to account for the different exposure lengths: near the maxima of 
both systems (phase 0.75 of J093010A and 0.63 of J093010B) 
the ratio is 44%, and near contact primary minimum (around 
phase 0.18 of J093010A and 0.0 of J093010B) it is 24%, show¬ 
ing considerable similarity across optical and near-infrared ob¬ 
servations. These consistent estimates provide a starting point 
for further modelling of the five stars. 


4. Modelling and Discussion 


To determine parameters for the eclipsing components of 
J093010, the two bi naries were modelled separately using the 
P HOEBE software (Prsa & Zwitteiil200^ . built upon the code 
of lWilson & Devinn^ ( 19711) . 

The starting assumptions for J093010A, based on its light 
curve, were that it is a detached binary in a circular orbit viewed 
nearly edge-on, with M\ > M 2 , T\ > T 2 , a third light con¬ 
tributing up to half the total flux, and a spot on one com po- 
nent. The mass ratio (q - 0.806, where iKoo et al.] (1201 4l) had 
found 0.841 by modelling) found from spectral disentangling 
was set as a non-adjustable parameter; the period was also fixed 
at 1.3055428 d (from the SuperWASP results) and the reference 
HJDq at 2456346.78443 (from iKoo et all ’s single primary min¬ 
imum timing); data was not binned or converted to phases prior 
to modelling. With the angle of inclination i set to 90°, the radial 
velocity curves were first fitted for a sin / (a proxy to the semi¬ 
major axis) and system velocity yo simultaneously, and the op¬ 
timum values found were not adjusted subsequently, since light 
curves are largely insensitive to these parameters. 

Pig. El shows the best fit obtain ed for the combination of 
our radial velocity data and that of iKoo et al.L Similar results 
were found when modelling our data sets separately: a sin i of 
5.765 + 0.014 Rq and 5.607 + 0.014 Rq for us and them respec¬ 
tively; and yo of -12.3 + 0.2 km s * and -10.4 + 0.2 km s *. 
Combining the data sets reduced the formal uncertainties on 
the two parameters in comparison with each data set considered 
alone (a sin / = 5.759+0.010 Rq and yo = -12.24+0.17 km s"'), 
and this combination of radial velocity data was used for the re¬ 
mainder of th e modellin g. 

To model 1^0 et aLf s light curves, a valid third light contri¬ 
bution needed to be i ncluded, m ore precise than the < 53% es¬ 
timated above. While l^o et al.l seem to have fitted their curves 
treating third light as a free parameter, we feared there might be 
strong correlations between this and other fitted parameters, and 



Phose 


Fig. 14. PHOEBE best model fit (red solid curve) to primary 
eclipse of J093010A in V, showing good match to eclipse width, 
shape and depth. 


so sought to constrain the value of the third light independently. 
Since the eclipses of J093010A appear total-annular, we could 
determine the phases of first, second, third and fourth contact 
((pi o ^ 4 ) quite precisely for the primary eclipse (see Pig. [TO. us- 
ing lKoo et al.l ’ s BV photometry, and then estimate the radii using 
these formulae from iHilditchl I2OOII) : 

7/? 

(02 - <^l) = (04 - <^ 3 ) = (1) 

Ina 

'IR 

(03 - 0l) = (04 - 02) = 2 ^- (2) 

On the assumption that the primary had the larger radius, 
this gave R\ - 0.832 + 0.018Rq and R 2 — 0.669 + 0.018Rq 
(where iKoo et afl found R\ - 0.757 + 0.008Rq and R 2 - 
0.743 + O.OlO RrT)). Then, we made use of an eclipse modelling 
guideline from lWilsonI (Il994l) : “the ratio of the depth of the an¬ 
nular eclipse to the light remaining in the total eclipse is approx¬ 
imately the square of the ratio of smaller to larger star radii”. We 
measured the depth of the annular (primary) eclipse in V and 
B (in PHOEBE flux units) from the light curves, and combined 
this with the ratio (f^)^ to estimate the flux which would remain 
in the total (secondary) eclipse, if there were no third light con¬ 
tribution: this would be the true contribution of the primary star 
(1.33 + 0.09 units in V and 0.56 + 0.04 units in B). The sec¬ 
ondary’s flux was given directly - and hence with greater pre¬ 
cision - as the depth of the total eclipse (0.275 + 0.011 units 
in V and 0.090 + 0.007 in B). By subtracting the primary and 
secondary fluxes from the maximum flux in the light curves, the 
tertiary’s contribution could then be estimated more precisely 
(0.90 + 0.09 units in V and 0.37 + 0.04 in B) as 36% + 4: some¬ 
what larger than the fitted values found bv iKoo et alJ of 31% in 
V and 32% in B. These values for third light were then fixed for 
the rest of the modelling. 

The Kopal potentials for the two eclipsing stars were 
then set to values which would reproduce the radii determined 
above (7.74 and 8.04 respectively), indicating strongly negative 
filling factors, as expected for a well-detached system. These 
values gave excellent matches to the widths of the eclipses 
(e.g. Eig. [T4l l. The effective temperatures of both components 
were also estimated from their B - V colours, using the indi¬ 
vidual fluxes found above: these were 4970^210^ for the pri¬ 
mary and 4470 + 160K for the better-determined secondary (by 
the same method, the third component was estimated to have 
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Fig, 15, Top: V band light curves for J093010A (data from 
iKoo et al.1 (120141) ') expressed in PHOEBE flux units, with best¬ 
fitting PHOEBE model curves overplotted (dashed red lines in¬ 
dicate model with no spot, solid blue lines indicate one-spot ex¬ 
ample model). Middle: residuals from no-spot model, showing 
out-of-eclipse semi-sinusoidal light variation unaccounted for. 
Bottom: residuals from one-spot model, largely correcting for 
these variations. 


Fig. 16. B band light curves, fits and residuals for no-spot and 
one-spot models for J093010A (see caption to Fig. dSl). 



r zi 


r 3 = 4900^3 jQ, consistent with the value found from its disen¬ 
tangled spectrum earlier). Using these values for T i and T 2 gave 
a model which had too shallow a primary eclipse and too deep 
a secondary eclipse; their ratio was therefore increased by in¬ 
crements until both eclipses were equally well reproduced: this 
occurred at Ti = 5185K and T 2 - 4325K. The angle of incli¬ 
nation i also required a slight decrease from 90° to match the 
depths of both eclipses simultaneously (hence a and were 
also fractionally adjusted). As a further check on these tem¬ 
peratures, we may note that the ratio of primary to secondary 
eclipse depths matched the correspondin g ratio of the stars’ sur¬ 
face brightnesses, as would be expected ( 'Wilsoniri994l) . 

As noted earlier, the apparently flat-bottomed primary 
eclipse could alternatively suggest R 2 > Ri, and so this possibil¬ 
ity was also explored in the modelling. Using Equations (1) and 
(2) with the roles of Ri and R 2 reversed, the components’ radii, 
third light contribution, Qi ,2 and temperatures were again esti¬ 
mated as described above; however, no appropriate ratio of tem¬ 
peratures could be found which would fit both eclipse depths un¬ 
less the uncertainties on ^ were exceeded. Moreover, the best 
fit produced an inconsistency between the eclipse depth ratio and 
the surface brightness ratio. In view of these problems, and the 
mass-radius mismatch (which would be challenging to explain 
in a well-detached system), we prefer the original model for the 
remainder of the analysis. However, improved photometry, es¬ 
pecially of the primary eclipse, would be desirable in future to 
confirm it more fully. 

The full optimal light curve fits and their residuals are shown 
in Figs.fTSlandfTbl It is apparent that, while the eclipses are very 
well-fitted by this model, the out-of-eclipse portions are poorly 


Fig. 17. Image of best-fit PHOEBE one-spot model for 
J093010A at phase 0.25. The large cool spot is visible on one 
limb of the (smaller) primary at quadrature. 


fitted (reduced= 17.1 in "V, 2.62 in B). This is almost certainly 
due to the presence of the assumed moving spot discussed ear¬ 
lier. To test this explanation, a single large cool spot was placed 
on the primary and its four parameters varied manually until a 
fair match to the light curve was achieved (colatitude 90°, lon¬ 
gitude 29°, radius 95°and temperature 0.9835 x mean primary 
temperature). The resulting light curve fits and residuals are also 
shown in Figs. [15] and [16] and an image of the model system in 
Fig- [13 These models reproduce the light variations much better 
(X^ = 8.86 in V, 1.18 in B). 

It should be emphasized that we do not regard this as the 
only possible spot configuration which would reproduce the ob¬ 
served light curve variations: a cool spot on the opposite side of 
the secondary also works passably well, and ho t spots on e ither 
component would also achieve similar results. iKoo et all used 
two cool spots with quite different parameters from ours, one 
on each component in their preferred model, to achieve a very 
close fit to the light curve; however, by appropriate placement of 
multiple spots, any light curve features whatsoever can be repro¬ 
duced, and we did not wish to include more components in our 
model than were justifiable. That said, it is very likely that there 
are multiple spots on this system, or non-circular spots, whose 
modelling would produce an even better fit to the light curve. 
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Table 2. System and stellar component parameters for J093010A 



Primary 

Secondary 

Orbital period (s) 

112798.90 + 0.16 

Semi-major axis (R©) 

5.762 + 0.010 

Mass ratio 

0.806 + 0.007 

System velocity (km s'*) 

-12.24 

+ 0.17 

Angle of inclination (°) 

88.2 

+ 0.3 

Kopal potential 

7.74 + 0.15 

8.04 + 0.18 

Filling factor 

-9.4 

-10.0 

Mass (Mq) 

0.837 + 0.008 

0.674 + 0.007 

Radius (Rq) 

0.832 + 0.018 

0.669 + 0.018 

Temperature (K) 

5185+25 

4325+2“ 

Surface gravity 

4.52 

4.62 

Bolometric luminosity 

5.66 + 0.02 

6.92 + 0.02 

Flux contribution (V) 

0.532 + 0.036 

0.110 + 0.004 

Flux contribution (B) 

0.552 + 0.040 

0.089 + 0.007 


but a method such as Doppler tomography would be required to 
determine their nature more rigorously. 

Table l^gives the final parameters determined for J093010A, 
using the no-spot model (there is almost no change to any of 
these under the one-spot model). The uncertainties on a, yo and 
i are formal errors generated by PHOEBE; the uncertainties on 
q, El 2 , Ri ,2 and Qi _2 were obtained from independent measure¬ 
ments as described above; and the uncertainties on M 12 were 
found using the formula prescribed in the PHOEBE manual: 


cr Ml 


= Ml 





cr M 2 - ^2 




\ 

q{q+ 1)/' 


(3) 

(4) 


J093010B was simpler to model. We started with the as¬ 
sumptions that it is a contact system (i.e. T\ - T 2 - 4700K, 
Qi = Q 2 ) in a circular orbit seen nearly edge-on, with Mi > M 2 , 
Ri > R 2 and no third light. The period was set to 0.2277138 d, 
the mass ratio to 0.397 + 0.006, and the reference HJDq to 
2456288.8780 (an average of the preferred values for radial ve¬ 
locity curves and light curves, since the continuous variation of 
J093010B’s light curves made estimation of the time of primary 
minimum less reliable from photometry alone). We first fitted 
the radial velocity curves alone for a sin / and yo, as we did for 
J093010A, and found the best fit shown in Eig. |7l which gave 
a sin / = 1.661 +0.009R£) and 70 = -11.3 + 0.7kms‘^. Here, we 
note that iKoo et alJ (1201 4ll also advanced a model for J093010B 
based purely on its light curve, which used a value for the mass 
ratio of 0.468 + 0.005 found via the “q-search” method: some 
distance from t he value fo und here. 

Then, using iKoo et alJ s BV light curves, and again adjusting 
a manually to keep a sin / at the value found earlier, / and Q were 
allowed to reach their optima simultaneously (the continuous 
variation of the light curve meant that the approach taken to find 
J093010A’s radii was impossible here). The best fits are shown 
in Eigs.[T8]and[T9l The Kopal potential corresponding to this fit 
(2.63+0.08) indicates a positive filling factor of 0.19 as expected 
for a contact system. The angle of inclination found is 86 ±4°, 
which we note is consistent with that found for J093010A, sug¬ 
gesting that the two binaries may share an orbital plane and 
perhaps originally fragmented out of a single protostellar disk. 
However, measurements of the longitude of the ascending node 
for both binaries would be required to prove a common plane. 

Again, a better fit can be achieved if a spot is modelled, to ac¬ 
count for the different heights of the maxima. An example spot¬ 



Fig. 18. V band light curves, fits and residuals for no-spot and 
one-spot models for J093010B (see caption to Eig. dsii. 



Fig. 19. B band light curves, fits and residuals for no-spot and 
one-spot models for J093010B (see caption to Eig.fTsTl. 


ted model is shown in Figs. [18] to |20j this uses a spot on the 
primary with colatitude 90°, longitude 60°, radius 110°and tem¬ 
perature 0.979 X mean primary temperature. As with J093010A, 
we emphasize that this is only one of several possible w ays to 
model the light variations using a single spot; iKoo et afl used a 
cool spot on the primary and a hot spot on the secondary instead. 
Our resulting parameters for the system, shown in Table |3| are 
not affected by the presence of the spot. The uncertainties on a, 
q, 70 . i and Mi^ were found as for J093010A; that for T was 
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Fig. 20. Image of best-fit PHOEBE one-spot model for 
J093010B at phase 0.25. The large cool spot is visible at the 
limb of the (larger) primary. 

Table 3. System and stellar component parameters for J093010B 



Primary 

Secondary 

Orbital period (s) 

19674.47 + 0.03 

Semi-major axis (R©) 

1.665 +0.012 

Mass ratio 

0.397 ± 0.006 

System velocity (km s"') 

-11.3 + 0.7 

Angle of incl. (°) 

86 + 4 

Kopal potential 

2.63 ± 0.08 

Filling factor 

0.17 

Mass (M©) 

0.86 + 0.02 0.341 + 0.011 

Radius (R©) 

0.79 + 0.04 

0.52 + 0.05 

Temperature (K) 

4700 + 50 

Surface gravity 

4.58 

4.53 

Bolometric luminosity 

6.20 + 0.05 

7.12 + 0.05 


estimated from the fitting to the disentangled primary spectrum; 
that for O is a formal output of PHOEBE; and those for /?i ,2 were 
found by setting a, q and Q to their extrema in appropriate com¬ 
binations, which may be expected to overestimate these errors 
somewhat. 

Eor both systems, dis tances we re estimated by a similar 
method to that outlined in iKoo et all absolute bolometric mag¬ 
nitudes (outputs from PHOEBE) were converted to abs olute V 
magni tudes using bolometric corrections as ta bulated in Elowed 
(Il996l) . and then combined using the formula (lHilditchll2001 ) 


Mv,,c,al - Mv ,2 = -2.5 logio(l + 


(5) 


Apparent V magnitudes of the maxim um brightness of the sys¬ 
tem were determined from iKoo et al.l ’s V light curves at their 
maxima for J093010B; for J093010A the contribution of the 
third star had to be removed first, using the component lig ht ra¬ 
tios obtained during analysis. The formula (lHilditchll200lll : 


V-Av - Mv,total ^ 5logiQd-5 (6) 

was then used to obtain distance estimates in parsecs (taking 
Ay - 0.03). Eor J093010A, d = 78 + 3 pc was found, and for 
J093010B , d - 11> ± 4 pc; the ranges overlap substantially at 
75-77 PC. IKoo et alJ found a consistent, if more uncertain value 
for J093010B of 77+9 pc, but a shorter distance of 66+7 pc for 
J093010A, probably due in the main to their handling of the third 
star’s contribution to system light. The similarity of our distance 
estimates lends further support to the reality of the association 
between the two binaries. 

An approximate distance was also estimated for the fifth star: 
its apparent V magnitude was found similarly to those for the 
binaries, from its fractional light contribution to J093010A and 



Fig. 21. Temperature-luminosity plot for primary of J093010A 
binary, compared with Dartmouth stellar evolution model 
isochrones. Each short curved line is a section of a 
model isochrone corresponding to masses between 0.829 and 
0.845 Mq i.e. the range of likely masses found from our mod¬ 
elling. The sections in the middle (black) have solar metallicity; 
those on the right (red) have [Ee/H]=0.2 and those on the left 
(blue) have [Ee/H]=-0.05. Ages of four isochrones consistent 
with the primary’s parameters are indicated. 


from the V light curves of iKoo et al.L Its absolute V magnitude 
was assumed to be that of a typical mid-main-sequence star with 
solar metallicity and T - 5000 + 100 K (combining the two 
independent effective temperature measurements obtained ear- 
lier), and was take n from Dartmouth stellar evolution models 
(iDotter et al.ll2008|P i The formula above then gave d - 81 pc 
for a 5 Gyr star, or c/ = 83+9 pc for 10 Gyr; both ranges are 
consistent with the 75-77 pc distance estimate found for the two 
binaries. Although less reliable than the binary distances, this 
finding adds some support to the claim that the fifth star is phys¬ 
ically associated with J093010. 

The parameters for J093010A are fully consistent with 
Dartmouth models. Plotting temperature against luminosity, 
each binary component overlaps with the range of values ex¬ 
pected for a star of its mass, assuming solar metallicity (Eigs.lTTI 
and l22l l. A very tiny negative value of [Ee/H] is p erhaps im- 
plied for each, but nothing like the -0.25 found by iKoo et all 
An optimal average age of 10.3 + 0.7 Gyr is indicated by this 
method. Then, plotting masses against radii, both components 
are consistent with solar metallicity isochrones between around 
8 and 10 Gyr (Eig. l23t : this implies an optimal average age of 
9.2+ 1.0 Gyr. 

Fig- EH also sho ws that J0930 10A’s binary components be¬ 
have very similarly to lTorres et al.l ’s collection of low-mass stars 
in deta ched eclipsing binaries with well-determined parameters 
(l2010h . The components of J093010B, however, do not sit on 
a single isochrone, presumably because of their interaction in 
contact: the primary has a smaller radius than its mass would 

^ http://stellar.dartmouth.edu/models/ 
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Fig. 22. Temperature-luminosity plot for secondary of J09301OA 
binary, compared with sections of model isochrones correspond¬ 
ing to masses between 0.667 and 0.681 Mq (see caption to 
Fig. |2T]for further explanation). 



Moss / solor mosses 


Fig. 23. Mass-radius diagram for binary components of J093010, 
compared with Dartmouth model isochrones for solar metallic- 
ity (dotted lines, corresponding to 1, 8, 9, 10 and 14 Gyr from 
bottom to top), and with empirical masses and radii for collected 
detached (red diamonds) and contact (blue triangles) eclipsing 
binary components. The binary components of J093010A are in¬ 
dicated by larger bold-face red diamonds, and those of J093010B 
by larger bold-face blue triangles. 


imply, assuming an age of 9 or 10 Gyr, while the secondary’s 
radius is much greater (and is inconsistent with any possible 
isochrone). Comparison with the para meters of short-period con - 
tact binary components collected by Ste^ien&Gazea^ (1201 2h . 
supplemented by two determined in lLohr et al.l ( 2014a ), suggests 
that this behaviour is entirely normal for systems of this config¬ 
uration; all the secondaries (those with masses below 0.7 Mq) 
lie in this region of the mass-radius diagram, with clearly in¬ 
flated radii, while the primaries in most cases have radii consis¬ 
tent with or even smaller than those expected from normal stellar 
evolution models. 


5. Conclusion 

After our initial discovery of J0 93010 as a probable doubly 
eclipsing quadruple system, and iKoo et al.l ' s follow-up obser¬ 
vations which confirmed J093010A as a double-lined spectro¬ 
scopic binary and revealed the presence of a possible fifth com¬ 
ponent, we have carried out further spectroscopic and photo¬ 
metric observations of the whole system. We have confirmed 
J093010B as a double-lined spectroscopic binary also, and con¬ 
firmed the association on the sky of the fifth component within 
J093010A. Radial velocities were measured for the first time for 
J093010B, allowing a more reliable determination of its mass ra¬ 
tio; the effective temperature of its primary, and of the assumed 
tertiary in J093010A, were also estimated from their disentan¬ 
gled individual spectra. Reanalysed SuperWASP photometry in¬ 
dicated the presence of a moving spot on one component of the 
binary in J093010A, with an apparent speed relative to its host 
star of approximately one rotation every two years. 

After modelling both eclipsing systems w e have obta ined pa¬ 
rameters for them which agree with those of iKoo et al.l in some 
areas, but disagree in several others, owing to different modelling 
approaches and available data. Specifically, we find a greater 
mass ratio for J093010B on the basis of our new spectroscopic 
data, and different radii for the binary components of J093010A 
as a result of our use of the total-annular eclipses in its light 
curves, which also allowed direct calculation of the third light 
contributed by the tertiary. Both binaries probably possess at 
least one large spot, though further observations would be re¬ 
quired to constrain any spot properties plausibly. 

On the basis of our temperature estimates from the disen¬ 
tangled component spectra, and from modelling, consistent dis¬ 
tances for both binaries of 73 and 78 pc are found. The fifth 
non-eclipsing star appears to have a temperature and spectrum 
similar to those of the primaries in the two eclipsing binaries; a 
consistent distance is also found for it if it is assumed to be a typ¬ 
ical main-sequence star. The system velocities of both binaries, 
and the average radial velocity of the fifth star, are all around 
-11 to -12 km s ’. These findings strongly support the claim 
that both binaries, and very probably all five stars are gravita¬ 
tionally bound in a single system. The consistent angles of in¬ 
clination of the two binary systems also imply likely formation 
by fragmentation from a single protostellar disk, without sub¬ 
sequent significant disruption to their orbital plane. The age of 
the system is around 9-10 Gyr, and it has approximately solar 
metal licity. 

This bright, close, highly unusual star system, containing a 
very short-period contact eclipsing binary, and a triple system 
including a low-mass detached eclipsing binary, would doubtless 
repay further investigation. 
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Table 1. Spectroscopic observations and heliocentric-corrected radial velocities for J093010. 


HJD 

Exp. time /s 

(A) Prim. RV /km s'* 

S Prim. RV 

(A) Sec. RV 

S Sec. RV 

(A) Ter. RV 

S Ter. RV 

(B) Prim. RV 

6 Prim. RV 

(B) Sec. RV 

S Sec. RV 

2456283.7679 

5 

87.5 

1.0 

-132.3 

1.2 

-11.8 

0.5 

28.9 

1.5 

-116.4 

1.5 

2456283.7704 

30 

90.1 

0.7 

-132.5 

1.1 

-11.1 

0.6 

40.4 

1.5 

-138.1 

1.5 

2456283.7710 

30 

88.6 

1.6 

-133.0 

1.3 

-11.1 

0.3 

39.2 

1.5 

-143.4 

1.5 

2456283.7716 

30 

89.7 

0.6 

-130.9 

0.5 

-11.3 

0.9 

44.5 

1.5 

-141.7 

1.5 

2456283.7749 

180 

88.3 

0.5 

-133.3 

0.9 

-11.4 

0.6 

51.5 

1.5 

-162.8 

1.5 

2456283.7773 

180 

88.7 

0.8 

-133.2 

0.3 

-11.0 

0.3 

52.6 

1.5 

-179.0 

1.5 

2456283.7796 

180 

88.4 

0.4 

-133.3 

0.3 

-10.98 

0.25 

55.9 

1.5 

-180.8 

1.5 

2456283.7853 

360 

89.1 

0.4 

-133.0 

0.8 

-11.2 

0.8 

69.2 

1.5 

-203.4 

1.5 

2456283.7897 

180 

89.0 

0.5 

-133.4 

1.2 

-11.5 

0.4 

71.1 

1.5 

-219.4 

1.5 

2456283.7921 

180 

88.7 

0.3 

-133.0 

0.9 

-10.58 

0.14 

74.1 

1.5 

-237.1 

1.5 

2456283.7979 

360 

86.8 

1.0 

-134.2 

2.0 

-13.0 

1.6 





2456292.7661 

60 

58.01 

0.18 

-102.3 

1.5 

-12.35 

0.18 

-64.6 

1.5 

124.7 

1.5 

2456292.7671 

60 

60.3 

0.7 

-101 

3 

-12.3 

0.6 

-63.3 

1.5 

125.4 

1.5 

2456292.7680 

60 

59.0 

0.6 

-99.6 

1.3 

-11.8 

0.5 

-66.7 

1.5 

132.3 

1.5 

2456292.7710 

10 

58.8 

1.6 

-98 

5 

-13.3 

1.1 

-79.4 

1.5 

147.3 

1.5 

2456292.7743 

10 

55.9 

1.3 

-102.5 

2.2 

-13.1 

0.3 

-84.1 

1.5 

163.5 

1.5 

2456292.7758 

60 

59.4 

0.8 

-106.4 

0.8 

-13.39 

0.19 

-88.6 

1.5 

170.3 

1.5 

2456292.7767 

60 

59.7 

0.8 

-106.1 

2.5 

-13.1 

0.5 

-88.7 

1.5 

173.2 

1.5 

2456292.7777 

60 

59.33 

0.12 

-107.54 

0.18 

-14.3 

0.7 

-93.9 

1.5 

176.0 

1.5 

2456292.7822 

360 

61.2 

0.5 

-106.6 

0.4 

-13.2 

0.4 

-100.4 

1.5 

197.7 

1.5 

2456292.7866 

360 

62.5 

0.6 

-107.6 

0.5 

-13.2 

0.5 

-107.0 

1.5 

213.8 

1.5 

2456292.7911 

360 

63.8 

0.4 

-110.5 

0.7 

-13.16 

0.18 

-109.8 

1.5 

229.5 

1.5 

2456292.7948 

60 

67.2 

1.5 

-113.0 

0.6 

-16.8 

1.8 

-118.7 

1.5 

237.3 

1.5 

2456292.7957 

60 

68.8 

2.3 

-110.2 

0.6 

-10.9 

0.7 

-116.4 

1.5 

239.8 

1.5 

2456292.7967 

60 

69.3 

0.8 

-112.5 

0.9 

-12.8 

0.9 

-114.4 

1.5 

246.6 

1.5 

2456292.7992 

300 

66.4 

0.7 

-112.4 

0.8 

-13.6 

0.4 

-116.0 

1.5 

247.1 

1.5 

2456292.8029 

300 

67.7 

0.4 

-112.2 

0.9 

-12.83 

0.13 

-119.0 

1.5 

247.5 

1.5 

2456292.8066 

300 

68.49 

0.06 

-114.9 

0.9 

-13.1 

0.3 

-117.5 

1.5 

247.6 

1.5 

2456294.7665 

60 

-91.7 

0.8 

96.0 

1.7 

-9.05 

0.13 

67.1 

1.5 

-209.4 

1.5 

2456294.7674 

60 

-91.4 

1.4 

96.3 

2.2 

-9.1 

0.3 

66.6 

1.5 

-196.4 

1.5 

2456294.7684 

60 

-93.6 

0.6 

94.4 

1.1 

-9.6 

0.5 

64.0 

1.5 

-197.0 

1.5 

2456294.7707 

120 

-94.4 

0.6 

94.9 

0.9 

-10.25 

0.10 

67.8 

1.5 

-178.5 

1.5 

2456294.7724 

120 

-94.4 

0.7 

95.9 

0.3 

-9.67 

0.15 

58.6 

1.5 

-179.1 

1.5 

2456294.7740 

120 

-95.4 

0.4 

95.7 

0.4 

-10.40 

0.06 

56.1 

1.5 

-171.7 

1.5 

2456294.7783 

110 

-95.59 

0.14 

97.0 

1.2 

-10.07 

0.21 

45.0 

1.5 

-143.4 

1.5 

2456294.7799 

no 

-97.0 

0.4 

96.4 

1.1 

-10.8 

0.5 

37.4 

1.5 

-143.4 

1.5 

2456294.7814 

no 

-96.6 

0.5 

96.6 

1.2 

-10.22 

0.20 

37.4 

1.5 

-126.2 

1.5 

2456294.7829 

no 

-97.8 

0.3 

96.7 

0.3 

-10.69 

0.09 

32.7 

1.5 

-128.3 

1.5 

2456294.7845 

no 

-97.2 

0.4 

98.1 

0.6 

-10.2 

0.3 

31.0 

1.5 

-108.3 

1.5 

2456294.7860 

no 

-98.2 

0.5 

96.9 

1.5 

-10.85 

0.08 

19.9 

1.5 

-95.7 

1.5 

2456294.7885 

no 

-99.0 

0.4 

96.5 

0.7 

-11.03 

0.12 

13.6 

1.5 

-76.9 

1.5 

2456294.7901 

no 

-99.2 

1.0 

97.3 

1.0 

-11.02 

0.24 

13.5 

1.5 

-73.1 

1.5 

2456294.7916 

no 

-98.8 

0.4 

97.8 

1.0 

-10.53 

0.21 

7.9 

1.5 

-56.5 

1.5 

2456294.7931 

no 

-100.02 

0.16 

96.9 

0.4 

-11.25 

0.24 

-1.4 

1.5 

-41.4 

1.5 

2456294.7947 

no 

-99.4 

0.4 

98.3 

0.8 

-10.61 

0.22 

-2.3 

1.5 

-35.0 

1.5 

2456294.7962 

no 

-100.3 

0.4 

97.7 

1.4 

-11.39 

0.20 

-5.2 

1.5 

-29.5 

1.5 

2456294.7986 

60 

-99.7 

0.6 

100.5 

1.3 

-10.5 

0.6 

-16.1 

1.5 

-3.7 

1.5 

2456294.7996 

60 

-100.5 

1.1 

98.7 

1.6 

-10.8 

0.3 

-10.2 

1.5 

-20.6 

1.5 

2456294.8005 

60 

-102.0 

0.7 

99.6 

0.5 

-11.9 

0.3 

-17.3 

1.5 

-5.8 

1.5 


M. E. Lohr et al.: The doubly eclipsing quintuple low-mass star system J093010 



